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In this work, we report on recent analysis of three-loop models of neutrino mass with dark matter. 

We discuss in detail the model of Krauss-Nasri-Trodden (KNT) P|, showing that it offers a viable 
solution to the neutrino mass and dark matter problems, and describe observable experimental 
signals predicted by the model. Furthermore, we show that the KNT model belongs to a larger class 
of three-loop models that can differ from the KNT approach in interesting ways. 


I. INTRODUCTION 

Models with radiative neutrino mass are of significant experimental interest. The inherent loop-suppression 
in such models allows the new physics responsible for neutrino mass to be lighter than otherwise expected. Such 
light new physics can be within experimental reach, either directly, through collider experiments, or indirectly, 
via e.g. searches for lepton flavor violating (LFV) effects. The loop suppression becomes more severe as the 
number of loops increases. Thus, models with three-loop masses are particularly interesting, as they generically 
require new physics at or around the TeV scale. 

He we present a class of models with radiative neutrino mass at the three-loop level [H-Q . We focus primarily 
on the KNT model [l| and report recent analysis showing that the model satisfies LFV constraints, such as 
/i —> e -|- 7 , and fits the neutrino oscillation data. Furthermore, the model contains a viable candidate for the 
dark matter (DM) in the universe, in the form of a light right-handed (RH) neutrino. We also show that a 
strongly first order electroweak phase transition can be achieved with a Higgs mass of ~ 125 GeV, as measured 
at the LHC ill- The model contains new charged scalars and we discuss their effect on the one-loop Higgs 
decay to neutral gauge bosons. Afterwards, we show that the KNT model belongs to a larger class of related 
three-loop models and briefly outline their features. 


II. THE KNT MODEL 

The KNT model [I| is an extension of the SM with three RH neutrinos, W, and two electrically charged 
scalars, and 5 ^, that are singlet under the SU{2)l gauge group. In addition, a discrete Z 2 symmetry is 
imposed on the model, under which {S' 2 , Ni} {—S 2 , —Ni}, and all other fields are even. This symmetry plays 
two key roles, preventing a tree-level coupling between Nn and the SM Higgs, which would otherwise induce 
tree-level neutrino masses, and ensuring that the lightest fermion Ni is a stable DM candidate. The Lagrangian 
reads 


C- = CsM + {fal3LlCiT2LgS^ + SiaNiS^laR + ^niNiN^Ni + h.c} -¥{^,81,82), 


( 1 ) 


where La is the left-handed lepton doublet, fa/s are Yukawa couplings which are antisymmetric in the generation 
indices a and /3, rriNi Majorana RH neutrino masses, C is the charge conjugation matrix, and V ($, S'!, S' 2 ) 

is the tree-level scalar potential, which is given by 


V(d>, Si, 2 ) = A ( 1 ^ 1 ")" - 1 ^ 1 " + m^S^Si + ml8;82 + Ai8*8i |$|" + A2S2*S2 

+ y (^i*-5i)' + y (S 2 *S 2 )' + r7i2Si*SiS2*S2 + {A,SiSiS 2 *S 2 * + h.c} . 


Id>l 


( 2 ) 


Here $ denotes the SM Higgs doublet. 
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FIG. 1: The three-loop diagram that generates the neutrino mass. 


A. Neutrino Mass and LFV 


The neutrino mass matrix elements, arising from the three-loop diagram in Fig. [TJ are given by 

Xsme^me^, 




ms- 


■faifpkgzjOkjF {m%jrri%^,rnljrnl^ , 


( 3 ) 


where p, k{= e, p, t) are the charged leptons flavor indices, i = 1, 2,3 denotes the three RH neutrinos, and the 
function F is a loop integral which is 0(1) It is interesting to note that, unlike the conventional seesaw 
mechanism, the radiatively generated neutrino masses are directly proportional to the charged lepton and RH 
neutrino masses, as well as being loop-suppressed. 

The Lagrangian © induces flavor violating processes such as la —>■ 'yip for > rnt,^, generated at one loop 
via the exchange of the charged scalars 3^2 ■ The branching ratio of such process is given by 


B{la 


3847r 


/.a/: 


k/3 


^Si 



( 4 ) 


with K ^ a,f3, aem being the fine structure constant and F 2 {x) = (1 — 6a; -I- 3a;^ -I- 2a;^ — 6a;^ lna;)/6(l — a;)^. 
For the case of la = Ip = P, this leads to a new contribution to the muon anomalous magnetic moment, 6a^. 


In our scan of the parameter space of the model we impose the experimental bound B (/r —>■ ey) < 
5.7 X 10“^^ [^; and use the allowed values for the neutrino mixing 5^2 = 0.320to;oi7i ^23 = 3.43to;o3j Sia = 
0.025)iQ oo|, and the mass squared difference = 2.55^0 09 x 10“^ eV^ and = 7.621°tjg x 10“°eV^ 

i- 


B. Dark matter 


An immediate implication of the Z 2 symmetry is that that lightest right handed neutrino, iVi, is stable, and 
hence a candidate for dark matter (DM). The Ni number density get depleted through the process NiNi —>• lalp 
via the t and u channels exchange of S^- In the non-relativistic limit, the total annihilation cross section is 
given by 


CkNiNi Vr 




48n {ml^ + 


( 5 ) 


with Vr is the relative velocity between the annihilation Ai’s. The relic density after the decoupling of Ni can 
be obtained by solving the Boltzmann equation, and it is approximately given by 




1.28 X IQ-^ / mjvi Y + 
Ea ./3 l5ia5j|';3p ^135 GeV/ l + 


( 6 ) 


where < ~ 6/a;/ ~ 6/25 is the thermal average of the relative velocity squared of a pair of two Ni particles, 

Mpi is planck mass; and g*{Tf) is the total number of effective massless degrees of freedom at the freeze-out 
temperature T/ ~ mArj/25. 
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FIG. 2: The charged scalar masses msi (red) and ms^ (green) versus the lightest RH neutrino mass, where the consistency 
with the neutrino data, LFV constraints and the DM relic density have been imposed. 


In Fig. [21 we plot the allowed mass range for , msj that gives the observed dark matter relic density 
Q. As seen in the figure, the neutrino experimental data, the bound on LFV process, combined with the relic 
density, seems to prefer for large regions of parameter space. However, the masses of both the DM 

and the charged scalar 5^ are bounded from above, with < 225 GeV and mgj < 245 GeV, respectively. 


C. Electroweak Phase Transition 

Although the SM has all the qualitative ingredients for electroweak baryogenesis, the amount of matter- 
antimatter asymmetry generated is too small. One of the reasons has to with the fact that the electroweak 
phase transition (EWPT) is not strongly first order, which is necessary to suppress the sphaleron processes 
in the broken phase. The strength of the EWPT can be improved if there are new scalar degrees of freedom 
around the electroweak scale coupled to the SM Higgs, which is the case in the KNT model. 

The investigation of the scalar effective potential reveals that, within the allowed parameter space of the 
model, the strength of the electroweak phase transition (EWPT) can be hrst order Q. We find that if the 
one-loop corrections to the Higgs mass are sizeable, then the condition v{Tc)/Tc > 1 for having a first order 
EWPT can be realized while keeping the Higgs mass around 125 GeV. The reason for this being that the extra 
charged singlets affect the dynamics of the SM scalar field VEV around the critical temperature M- In Fig. 
|3l we show the plot for v{Tc)/Tc versus the critical temperature and one observes that a strongly first order 
EWPT is possible while the critical temperature lies around 100 GeV. 



Tc (GeV) 


FIG. 3: The critical temperature is presented versus the quantity VcjTc. 
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Ecm {GeV) 

Selection cuts 

250 

70 GeV < El <11Q GeV, 70 GeV < M,,,, < 220 GeV, 
Mmiss < 120 GeV, 0.5 < 7?e < 2, — 2 < ry,, < —0.5. 

350 

90 GeV <Ei< 165 GeV, 100 GeV < Me,,, < 280 GeV, 
Mmias < 200 GeV, 0.5 < pe < 3, — 2.5 < r;,, < 0. 

500 

120 GeV <Ei< 240 GeV, 300 GeV < Me,„ < 480 GeV, 
Mmiss < 300 GeV, 0.5 < pe < 3, — 3 < Pf^ < 0. 

1000 

El < 70 GeV, Me.„ < 140 GeV, Mmiss > 750 GeV 

, 0.1 < pe < 0.8, — 0.8 < p,j, < —0.1. 


TABLE I: Relevant cuts for the process e"’"e —>■ Emiss +e pT for different CM energies. Here Ef, and ge are the charged 

lepton energy in its pseudo-rapidity, is the electron-muon invariant mass and M^iss is the missing invariant mass. 


D. KNT at High Energy Lepton colliders 

Since the RH neutrinos couple to the charged leptons, one excepts them to be produced at e“e+ colliders, 
such as the ILC and CLIC, with a collision energy ,/s of few hundreds of GeV up to a TeV. If the produced 
pairs are the form A^ 2 , 3 -^ 2,3 or A^iA^ 2 , 3 , then N 2,3 will decay into a charged lepton and 5^, and subsequently 
5"^ decays into Ni and a charged lepton. In addition, the SM neutrinos will be also produced via the decay of 
the charged scalar 51. Here we concentrate on he process e“e+ —>■ + Emiss [HI, where the missing energy 

corresponds to any state in the set Smiss C Ve^T, NiNk] i,k = 1, 2,3}. The total 

expected cross section of the processes e“e+ —>• e“/i+ + Emiss is represented by , while cr{£miss) denotes 
the cross section of different subprocesses. As a benchmark, we consider 

= -(4.97 + H.41) X 10"^ fer = 0.106 + i0.0859, ^ = (3-04 - 14.72) x 10"®, 

/ 0.2249+ i0.3252 0.0053 + 10.7789 0.4709+ zl.47 \ 

g^„ = 10-2x 1.099+ a.511 -1.365 -il.003 0.6532-10.1845 , 

\ 122.1 + 1178.4 -0.6398- 10.6656 -10.56 + 168.56 J 

mN, = {162.2 GeV, 182.1 GeV, 209.8 GeV}, ms, = (914.2 GeV, 239.7 GeV}, (7) 

We use LanHep and CalcHep to simulate our model and generate the differential cross section and the relevant 
kinematic variables for different CM energy: Ecm =250 , 350, 500 GeV andl TeV, with unpolarized beams at 
first; and then we consider the polarized beams with P (e“, e’*') = [—0.8, +0.3] and/or P (e~, e’*') = [+0.8, —0.3]. 
Imposing the appropriate cuts given in Table-jT] we show in Fig. 0] the dependance of the significance on the 
accumulated luminosity with and without polarized beams for the considered CM energies. We clearly see that 
for a polarized beam, the signal can be observed even with relatively low integrated luminosity. For example, at 
Ecm = 250 GeV, the 5 tr required luminosity is 150 fb~^ for polarized beam as compared to 700 fb~^ without 
polarization. 


III. A CLASS OF THREE-LOOP MODELS WITH DARK MATTER 

The KNT model is a simple theory that is capable of simultaneously addressing the neutrino mass and DM 
problems. It predicts new physics that can be probed at colliders and precision experiments. Within the KNT 
model, the DM is part of a Z 2 -odd sector that propagates in the inner loop of the neutrino mass diagram (see 
Figured]). One might wonder whether generalizations of the KNT model exist. It is well known that the seesaw 
mechanism can be generalized to a type-III variant with SU{2) triplet fermions [l^ . and a further variant with 
quintuplet fermions [l3l |. It is perhaps not surprising, therefore, to learn that the KNT model forms part of a 
larger class of theories that similarly achieve neutrino mass at the three-loop level, with DM propagating in the 
inner loop. 

The generalized models require that the SM be extended to include the charged singlet scalar ^ (1,1, 2), 
a new scalar S 2 , and three real fermions Np. However, instead of taking S 2 and Np as SU{2)p singlets, the 
variant models employ larger multiplets. The resulting models share key similarities with the KNT model, with 
the neutrino mass diagram retaining the exact form in Figured] and the DM remaining as the lightest fermion 
Np. However, there are also some interesting differences, relative to the KNT model. Firstly, the DM is now 
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FIG. 4: The significance versus luminosity at different CM energies within the cuts defined in Tableflji with (left) and 
without (right) polarized beams. The two horizontal dashed lines represent S — Z and <S = 5, respectively. 


part a larger multiplet with non-trivial gauge interactions. This modihes expectations for the DM mass and 
allows new DM interactions. Secondly, there are interesting consequences for the Z 2 symmetry, which depend 
on the specific details of the model, as we now outline: 

• The triplet model: Here S 2 ~ (1,3,2) and Nr ~ (1,3,0) are taken as SU{2)l triplets [l^. Similar to the 
KNT model, the triplet model requires a Z 2 symmetry with action {S 2 , Nr} —>■ {—/S' 2 , —Nr} to prevent 
tree-level neutrino masses and ensure a stable DM candidate. The DM is the neutral component of the 
lightest triplet fermion and has extra Yukawa interactions mediated by gia 7 ^ 0 in Eq. ((TJ. Consequently 
its mass increases to ~ 3 TeV [l3 |. 

• The quintuplet model: This model employs the multiplets S '2 ~ (1,5,2) and Nr ~ (1,5,0) [1^. Unlike 
the KNT model and the triplet model, the quintuplet variant does not require a Z 2 symmetry to prevent 
tree-level neutrino mass. It is therefore a viable model of radiative neutrino mass, independent of DM 
considerations. Interestingly, the Z 2 symmetry {S 2 ,iV;j} —>• {—S 2 ,—Nr} is broken by a single coupling 
A in the model; for technically-natural values of A <C 1, one obtains a long-lived DM candidate, while in 
the limit A —7 0 the Z 2 symmetry becomes exact and the DM is absolutely stable [l^. The quintuplet 
model is therefore a model of radiative neutrino mass, with or without DM, depending on the region of 
parameter space considered. Due to the gauge interactions, the case with DM requires ~ 10 TeV. 

• The septuple! model: This case has S 2 (1,7,2) and Nr ~ (1,7,0), and is of particular interest [l^ . 
Similar to the quintuplet model, the Z 2 symmetry is not required to prevent tree-level neutrino mass. 
However, different from the quintuplet model, the septuple! model automatically possesses an accidental 
Z 2 symmetry with action {S' 2 , iVi^} — >• {— S 2 , —Nr}, and therefore always contains a stable DM candidate. 
This gives a common description for neutrino mass and DM without requiring any new symmetries. The 
DM is relatively heavy, with ~ 20 — 25 TeV, due to the Sommerfeld enhancement induced by SU{2) 
gauge boson exchange [l^ . 

Detailed studies show that the variant models have large regions of viable parameter space that fit the neutrino 
oscillation data, reproduce the observed DM relic density, and satisfy experimental constraints from, e.g, LEV 
searches [idl - H^ . In each of the variant models, the DM is a heavy Majorana fermion and the mass for S 2 
should obey M 2 > M-^m, placing both multiplets beyond the reach of current colliders. However, the scalar Si 
can retain an 0(100) GeV mass in all cases, offering the best chance for testing the models at colliders. The 
models with larger SU (2 )r multiplets also generate sizable contributions to LEV processes like /i —>■ e -I- 7 . In 
fact, these are typically enhanced relative to the KNT case, improving the prospects for testing the model in 
future LEV searches [lil - H^ . Due to the new gauge interactions for the DM, direct-detection prospects also 
improve for the models with larger multiplets |14l4l6l |. 

There are additional models of neutrino mass with DM that are related to the KNT model [I3- For ex¬ 
ample, one can replace the internal SM leptons in Figure [T] with down-type quarks, cr^r —7 cIr^r, and take 
Si ~ (3,1,2/3) and S 2 ~ (3,1,—2/3) as colored multiplets. Further colored variants are also possible m. 
Alternatively, one may replace the real (Majorana) fermion Nr with a complex (Dirac) fermion and consider 
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FIG. 5: Three-loop diagram for neutrino mass in a variant model with Si ~ (1, 2,1), S 2 (1, 2, 3) and F ~ (1, 2, —1). 
The beyond-SM multiplets Si ,2 and F are taken odd under a Z 2 symmetry and the dark matter belongs to the inert 
doublet Si. 


a modified topology for the loop-diagram, with the DM being an inert scalar (see Figure [5] for an example). A 
systematic description of these variant models appears in Ref. da. 


IV. CONCLUSION 

Models with radiative neutrino mass offer a promising way to experimentally discern the new physics re¬ 
sponsible for the origin of neutrino mass. Three-loop models are particularly interesting as the new physics is 
generically expected at or around the TeV scale. Furthermore, such models can provide viable DM candidates, 
thereby solving both the neutrino mass and DM problems. Here we reported recent analysis of the KNT model, 
which show the model to be a viable theory for neutrino mass and DM that can be tested at collider experiments. 
We also showed that the KNT model belongs to a larger class of three-loop models that can generate similarly 
interesting observable effects. 
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